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An algorithm for controlling the tempering process for automobile glass is simulated. It is shown that the me¬ 
chanical properties of the glass produced can still be improved. The effectiveness of using modeling for deve¬ 
loping corrective actions in the production of tempered glass is validated. 
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Statistical analysis of tempered glass has revealed a sub¬ 
stantial variance of the mechanical properties. The variation 
of the character of breakage reaches 21 - 30%. These charac¬ 
teristics largely depend on the tempering process, as is con¬ 
firmed by regression models describing the dependence of 
the indicators on the tempering regime [1]: 
maximum number fragments: 

y x = 442.3 - 26.6x 13 - 316.4x lg - 8.0x 19 ; (1) 

largest length of fragments (mm): 

y 2 = - 196.5 + 6.0x 13 + 154.8x lg + 1.6x 19 + 0.25x 6 ; (2) 

minimum number of fragments: 

y 3 = 140.4 - 0.4x u - 102.4x lg - 2.5x 19 + 0.4x 9 , (3) 

where x 13 is the number of flows; x 18 ) interval 2 left-hand; 
x 19 ) height of punch, mm; x 6 ) vault temperature in chamber 
2 in zone 2 of the tempering furnace, °C; x n ) temperature in 
chamber 4 in zone 12, °C; x 9 ) temperature in chamber 3 in 
zone 11, °C. 

The quality of the glass produced can be increased by 
controlling the tempering regime. Control is a multicrite- 
rional decision-making problem. A decision made for cor¬ 
recting a regime is evaluated by a set of criteria characteriz¬ 
ing the mechanical properties of the glass, as determined by 
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the requirements of GOST 5727-88. The problem of control¬ 
ling the technological tempering process can be formulated 
as follows: 

minimize the linear form of the criterion 

K = Cj jq + C 2 y 2 + C 3 y 3 , (4) 

where C { and C 2 are coefficients, under the constraints 

Ti^Ti, pr ; (5) 

y 2 — Jh, pr ’ ( 6 ) 

Jb^pr- 

All regime variables x t must fall within admissible limits: 

*«, min max' ( 8 ) 

The criterion (4) is determined by the requirements im¬ 
posed on the character of the breakage of the glass during the 
tests. The maximum number of fragments y l and their largest 
length y 2 must not exceed prescribed values >q pr and jApr- 
The inequality (7) limits the minimum number of fragments 
during tests which must not be less than y 3 pr . 

The inequality (8) reflects the need to pick a regime x, 
from the region prescribed by the system of limitations. The 
regime variables x, will be found in the region described by 
the regression equations (1) - (3). 

The problem (4) - (8) formulated above is a typical line¬ 
ar-programming problem and can be solved by well-known 
methods. 
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Fig. 1. Form of the penalty functions. 


In solving the problem of controlling the tempering pro¬ 
cess the type of article being produced, which influences the 
value of the factor variable x u (number of flows), is consi¬ 
dered to be given. During the production process, perturba¬ 
tions can change the coefficients in the regression equations, 
which can make the constraints (5) - (8) incompatible, and 
the problem formulated will not have a solution. 

In this connection a control algorithm where the “rigid” 
constraints (5) - (7) are replaced by “non-rigid” constraints 
by reducing the mathematical programming problem to solv¬ 
ing an unconditional minimization problem using the method 
of penalty functions [2], The penalty function F of the fol¬ 
lowing form is chosen as the control criterion: 

F=\ abs (max (j-j - y hpr , 0)) —-— + 

y 1, max 

abs (min (v 3 -y 3 , 0))--— + 

y 3, min 

X 2 abs (max (v 2 - y 2 pr , 0)) --, (9) 

T2, max 

where X are weighting factors for the terms of the penalty 
functions and set equal to 0.33. 

The terms in the penalty function (9) are normalized by 
dividing each one by the admissible value of the indicator 
(GOST 5727-88), which made it possible to write the pe¬ 
nalty in the form of an additive function of dimensionless 
quantities with the coefficients X. The tenns of the penalty 
function are piece-wise smooth functions, shown in Fig. 1. 

The first and third terms in the expression (9) have the 
form of the penalty function A. The penalty is imposed when 
the indictor y exceeds a prescribed value v pr . The maximum 


TABLE 1. 


Value 

rarameier 

coded 

prescribed 

admissible 

Number of fragments: 




maximum 

y 1 

235 

400 

minimum 

k3 

102 

40 

Largest length of fragments 

k2 

39 

75 


number of fragments and their largest length are such indica¬ 
tors. The second term in the expression (9), limiting the mi¬ 
nimum number of fragments during the tests, has the form of 
the function B. A penalty is imposed in cases when the num¬ 
ber of fragments becomes less than a prescribed amount 
y 2 pr . The slope of the functions determines the rate of impo¬ 
sition of the penalty; it is established by adjusting the coeffi¬ 
cients X. Figure 1 shows an example of a change of the value 
of the penalty in proportion to the deviation of the parame¬ 
ters from the prescribed value (y — y pt ). 

The control problem consists in searching for a temper¬ 
ing regime where the penalty functions (9) assumes the mini¬ 
mum value. The penalty function (9) equals zero when each 
term is zero. The search for the optimal regime was made 
using a coordinate descent method [2], The region of search 
with respect to each regime variable (coordinate) was estab¬ 
lished by setting their range of variation. This range was 
taken to be the range of variation of the regime variables un¬ 
der industrial process conditions. The size of a stepped 
change of the regime variables during the search was 
adopted taking into account the accuracy of their measure¬ 
ment in a manner so that several steps could be taken in a 
fixed range of variation of the regimes. 

The algorithm described above was used to model the 
control algorithm used for tempering glass. The experiment 
was performed by using retrospective data on the operation 
of the tempering furnace during the production of colorless, 
curved, side window glass for automobiles. The objective of 
the experiment was to evaluate the possibility of increasing 
the quality of the glass produced by optimizing the temper¬ 
ing regime. 

In the modeling, requirements were prescribed for the 
quality of the tempered glass. These requirements are given 
in Table 1. The search range for the regime variables was 
taken as their range of variation under the conditions of com¬ 
mercial operation of the tempering furnace. 

The stepped change of the regime parameters in the 
search was chosen to be 10% of the range of variation. The 
range of variation of the regime variables in the modeling 
process are presented in Table 2. 

The modeling confirmed that the quality of the glass pro¬ 
duced can be increased. The comparative data of the sirnula- 


TABLE 2. 




Value 


Stepped 

Regime (factor) variable 

coded 

minimum 

maximum 

change 

Vault temperature in cham¬ 
ber 2 in zone 2 

Bottom temperature: 

x 6 

599 

645 

5 

in chamber 3 in zone 11 

Xg 

599 

654 

5 

in chamber 4 in zone 12 

x u 

599 

656 

5 

Interval 2 left-hand side 

X 1S 

0.35 

0.4 

0.00025 

Punch height 

X 19 

5 

18 

0.065 
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Experimental data 

Fig. 2. Number and size of the glass fragments in tests: 1 ) maxi¬ 
mum number of fragments (actual number); 2 ) maximum number of 
fragments (results of modeling); 3 ) minimum number of fragments 
(actual number); 4 ) minimum number of fragments (results of mod¬ 
eling); 5 ) maximum length of fragments (actual length); 6 ) maxi¬ 
mum length of fragments (results of modeling). 


tion of the control algorithm with manual insertion of the 
tempering process are presented in Table 3. As one can see 
from Table 3 the proposed control algorithm makes it possi¬ 
ble to produce tempered glass with prescribed quality re¬ 
quirements (see Table 1). 

The results of the tests of the glass for breakage are pre¬ 
sented in Fig. 2. Evidently, the maximum and minimum 
number of fragments and their largest length stabilize as 
compared with the manual insertion of the tempering process. 

The characteristics of the glass stabilize as a result of ad¬ 
justing the tempering regimes and holding them at the com¬ 
puted levels. Comparative data obtained on the tempering re¬ 
gimes by modeling with manual insertion of the process are 
presented in Table 4. 

The thermal regime of the tempering furnace varies very 
little, by no more than 6%. The punch height is substantially 
corrected to lower values. It is corrected by changing the 
number of glass process flows from 1 to 2 with the height set 
from 8 to 6 mm. The number of glass flows v 13 depends on 
the type of article produced. This factor was studied in the 
control problem as a planned prescribed quantity. 


TABLE 3. 


Parameter 


Coded 

value 


Modeling results 


Manual insertion 
of the process 


average standard average standard 
value deviation value deviation 


Number of frag¬ 
ments: 


maximum 

y\ 

230 

5.6 

161 

41 

minimum 


97 

2.5 

61 

9 

Largest length 
of fragments 


39 

0.2 

56 

10 


TABLE 4. 



Coded 

value 

Modeling results 

Manual insertion 
of the process 

Regime (factor) variable 

average 

value 

standard 

devia¬ 

tion 

average 

value 

standard 

devia¬ 

tion 

Vault temperature in 
chamber 2 in zone 2 

*6 

618 

5 

627 

11.5 

Bottom temperature: 






in chamber 3 in zone 11 

X 9 

649 

0 

612 

16.3 

in chamber 4 in zone 12 

*11 

600 

0 

614 

18.3 

Interval 2 left-hand side 

*18 

0.375 

0 

0.385 

0.02 

Punch height 

*19 

7 

1 

14.9 
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The simulation modeling performed showed that further 
improvement of the mechanical properties of the glass pro¬ 
duced is possible. It was also shown that modeling is helpful 
for developing corrective actions in quality management sys¬ 
tems in the production of tempered glass. 
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